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ABSTRACT 

We present the general relationship between phase correlations and the hierarchy of 
polyspectra in the Fourier space, and the new theoretical understanding of the phase 
information is provided. Phase correlations are related to the polyspectra only through 
the non-uniform distributions of the phase sum 6^ + • • • + 6k N with closed wave vectors, 
k± + • • • + Hn = 0. The exact relationship is given by the infinite series, which one can 
truncate in a consistent manner. The method to calculate the series to arbitrary order 
is explained, and the explicit expression of the first-order approximation is given. A 
numerical demonstration proves that the distribution of the phase sum is a robust 
estimator and provides an alternative statistic to search for the non-Gaussianity. 

Subject headings: cosmology: theory — large-scale structure of universe — methods: 
statistical 

1. Introduction 

Quantifying the cosmic fields, such as the density fields, velocity fields, gravitational lensing 
fields, temperature fluctuations in the cosmic microwave background, etc. is undoubtedly crucial 
to study the origin and dynamics of the structure in the universe. The structure of these fields are 
believed to be emerged from primordial random Gaussian perturbations, as most of the inflationary 
models naturally predict nearly scale-invariant Gaussian fluctuations (Guth & Pi 1982; Starobinskii 
1982; Hawking 1982; Bardeen, Steinhardt & Turner 1983). 

Even if the primordial perturbations are random Gaussian, the gravitationally nonlinear evo- 
lution produces non-Gaussianity in the cosmic fields. Quantifying the non-Gaussianity is not trivial 
since it depends on the full hierarchy of the higher-order correlation functions in real space, or of 
the polyspectra in Fourier space. First several members of such hierarchy can be observationally 
determined, with which only partial information on non-Gaussianity is quantified. Therefore, al- 
ternative statistics, such as the void probability function (White 1979), the genus statistic (Gott, 
Dickinson, & Melott 1986), the Minkowski functionals (Minkowski 1903; Mecke, Buchert &: Wag- 
ner 1994; Schmalzing & Buchert 1997), etc. which contain the information on the full hierarchy of 
higher-order statistics should be useful. 
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Non-Gaussianity is frequently termed "phase correlations" . This term reflects the fact that the 
Fourier phases of a random Gaussian field are randomly distributed without any correlation among 
different modes. Therefore, phase correlations, if any, obviously characterize the non-Gaussianity. 
However, what kind of phase correlations arise in a given non-Gaussian field have been far from 
obvious. Investigations along this line are quite limited in the literatures despite its importance, 
apparently because of the lack of theoretical guidelines. Most of the earlier work (Ryden & Gramann 
1991; Soda & Suto 1992; Jain & Bertschinger 1998) only assess the nonlinear evolution of phases 
in individual Fourier modes without statistics. Phenomenological studies of A-body simulations 
have revealed that the one-point phase distribution remains uniform even in non-Gaussian fields 
(Suginohara & Suto 1991), and that the phase difference between neighboring Fourier modes is 
non-uniformly distributed (Scherrer, Melott, & Shandarin 1991; Coles & Chiang 2000; Chiang 
2001; Chiang, Naselsky & Coles 2002; Watts, Coles, & Melott 2003). However, the meaning of the 
discovered phase correlations is obscure in those literatures. 

Since the hierarchy of the higher-order statistics contains statistically all information on the 
distribution (Bertschinger 1992), there should be some connection between phase correlations and 
polyspectra, which is the key to theoretically understand the phase correlations. Examining a toy 
model, Watts & Coles (2003) realized the importance of the phase sums with closed wave vectors in 
this connection, although they have never derived the exact relations. In this Letter, the connection 
in the general form is discovered for the first time. As a result, we will have much better theoretical 
understanding of the phase information than before. 



Although the real part Re/fc and the imaginary part Im/fc of the Fourier transform fk of a 
random field / are naturally the independent variables, one can also take their linear combinations 
fk = Re/fc + ilm/fc and ft = Re/fc — H m fk as another set of mutually independent variables. For 
calculational advantages, we use the latter choice. In this Letter, the reality of the random field / 
is assumed since most of the cosmic fields are real, although one can readily generalize the following 
analysis to the complex fields. Because of the reality condition, = /_&, f^s are actually not 
independent variables, and f^s of all modes k are taken as independent variables. 

It is useful to define the normalized quantity = fk/y/P{k), where P(k) = (\fk\ 2 ) is the 
power spectrum of the random field. The key technique to derive the relation between phase 
correlations and polyspectra is given by previous work (Matsubara 1995, 2003): the joint probability 
function P({ak}) of having particular set of is formally represented by 



2. Phase correlations and polyspectra 




(1) 



where (• • -) c indicates the cumulants, and "PG({ a fc}) is the multivariate Gaussian distribution func- 
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tion of variables {a k }. In the present case, ^({"fc}) exp (— \ Ylk a k a -k) , since the covariance 
matrix is (a k a k i) = <5j^ fe/ , where the symbol 5^ is defined by 5^ = 1 for k = and = for 
/ 0. The periodic boundary condition with boxsize V = L 3 is assumed. 

Since the polyspectra P( N \k±, . . . , fejv) are defined from the cumulants by 

</ fel ■ ■ ■ /* w >c = ^ 1 - JV/2 ^ 1 + ... +fcjvJ P W (^l, ■ ■ ■ , (2) 

the above formula (1) provides the relation between polyspectra and joint distribution of the Fourier 
coefficients. Expanding the exponential in equation (1), each term in this expansion consists of the 
products of polyspectra times derivatives of Vq. The derivatives of Vq are given by a simple 
polynomial of a k s times Vq. The general term in the expansion has the form 

£ Sl +k2+ ...5\ +k , + ... ■ ■ ■ P {N) (ki, k 2 , . . .) P (N) (k[,k' 2 , ...)••• H klk2 ... k , k , 2 ...V G , (3) 

k's 

with appropriate coefficients, where 

(AO/, , , \ P( N \k 1 ,k 2 ,...,k N - 1 ) 

p y >(ki, k 2 , ■ ■ ■ , fcjy-i) = — f — > (4) 

y/V^P^P^) • • • P(k N ^)P Wi + ■■■ + fciv-i|) 

are the dimensionless, normalized polyspectra of a k , and 



is a generalization of Hermite polynomials and is given by polynomials of a k s and 5^ 's. For 
example, 

H kl k 2 k 3 = a-fcia-fc 2 «-fc 3 ~ 5 fci+fc 2 a -fc3 - ^fca+fca^-fci ~ ^fc 3 +fci a -fc2> ( 6 ) 

and so on. Thus, the joint probability function "P({«fc}) is represented by Vq times infinite sum of 
products by a k s, S k 's, and normalized polyspectra. 

Next step is to transform the complex variable a k into the modulus \a k \ and the phase 9 k by 
Ofc = Wk\& ldk • This should be carefully done, since a k is considered independent on a* k . At this 
point, we restrict the wavenumber k in the upper half sphere (uhs), k z > 0, and the degrees of 
freedom in the lower half sphere is relabeled by the reality relation, a k = a*_ k for k z < 0. The 
mode k = is excluded which ensures zero mean of the original field /. The term (3) is accordingly 
relabeled, resulting in the sum of the products of a k s, otVs, S k 's, where k 6 uhs, and normalized 
polyspectra. With the above procedures, one can express the ratio 'P({afc})/'PG({ Q; fc}) in terms of 
the normalized polyspectra, the modulus \a k \ and the phase 6 k . The Jacobian of the transform 
from a k to (\a k \,6 k ) is the same for the probability functions V({a k }) and Vo({a k }). Therefore 
this meets our ends to relate the phase correlations and polyspectra, which is a completely new 
result. 



Practically, one needs to truncate the infinite series by a consistent manner. Fortunately, the 
non-Gaussianity generated by gravitationally nonlinear evolution is known to approximately follow 
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the hierarchical model of the higher-order correlations, in which the polyspectra have the 

order, p( JV ) ~ C[P(/c) Ar_1 ] (e.g., Bernardeau, Colombi, Gaztanaga & Scoccimarro 2002). This 
means ~ 0(e N ~ 2 ), where e ~ -y/ P(k) /V. Therefore one can evaluate the phase correlations in 
perturbative manner as long as the expansion parameter e is small. It is straightforward to perform 
the above procedure to express "P({|afc|, 9k}) in terms of normalized polyspectra to arbitrary order 
in e. In the lowest order approximation, only the normalized bispectrum p^ gives the term of 
order 0(e l ). The result is 

V({\a k \,6 k }) \{ d\a k \d6 k 

feguhs 

1+ \ a k 1 \\a k2 \\a kl+k2 \cos(e kl +e k2 -e kl+k2 )p^(k 1 ,k 2 ) 

fci,fc2Guhs 

d6 k 
27 



x J] 2\a k \e-\^ 2 d\a k \^. (7) 



feg uhs 

Higher-order terms can be similarly calculated, although they are somehow tedious. For example, 
in the second-order approximation, C(e 2 ), there appears the square of the first-order term, and 
terms like 

|afcil|afc2ll a fc3ll a fci+fc2±fc 3 l cos (^ fc i + ^ fc 2 ^ ^ fc 3 ~~ ^fci+fc 2 ±fc 3 ) (8) 
with appropriate normalized trispectrum or the product of normalized bispectra multiplied, and 
other terms which do not depend on phases. 

The phases always contribute to the probability distribution by the combination of the form, 
cos(6*fc 1 + • • • + 6k N ), with closed wavevectors: fci + • • • + fcjv = 0. This is generally true because the 
phase dependence in equation (3) is the exponential of the sum of phases, and the probability is 
the real number so that taking real parts gives the cosine function. The reason that phase correla- 
tions exist only among modes with closed wavevectors comes from the translational invariance. In 
equations (7) and (8), wavenumbers are restricted to the uhs so that the modes in the lower half 
sphere are relabeled by 9 k = —6- k . 

The moduli |afc|'s are easily integrated in the first-order approximation of equation (7), result- 
ing in 

/ — uhs i , — \ 3 uhs 

V{{9 k }) oc 1 + ^ J> (3) (fc, k) cos(2# fc - 9 2k ) + ( ^ ) E P (3) ( fc ' k ') cos (^ + 6 U ~ ( 9 ) 

k ^ ' k+k' 

The practically useful relations between phase correlations and the bispectrum are obtained by 
further integrating some phases in equation (9). One obtains 

r(e k , e 2k ) oc i + ^V 3) (fc, k) cos(20 fc - e 2k ), (10) 

3/2 

r(0k, o k ,,e k+k ,) oc i + ^ (3) (fc, k') cos (9 k + e w - o k+k ,), (ii) 
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where k / k' '. These are the explicit forms of the relation between phase correlations and the 
bispectrum in the first-order approximation. We find that the distribution of the 'phase sum' 
9 k + 9 k i — 9 k+k > is determined only by the normalized bispectrum at the first-order level, although 
higher-order normalized polyspectra can contribute in general. The higher-order calculations show 
that the distribution of the phase sum 9 kl + ■ ■ ■ + 9 kN for the modes with closed wavevectors 
k\ + • • • + &7v = is determined by normalized polyspectra of order 3 to N in the lowest-order 
approximation, where the identification 9 k = —9^ k is understood. It was vaguely suggested that 
there is some relationship between the phase sum and polyspectra based on a particular non- 
Gaussian model by Watts & Coles (2003). We now find the explicit relationship between them in 
general non-Gaussian fields. 

If we further integrate all phases but one particular 6 k , the one-point probability function of a 
phase is uniform, V(9 k ) = 1/27T, which is consistent with the previous TV-body analysis (Suginohara 
& Suto 1991). This conclusion does not depend on the first-order approximation, since a single 
wavevector can not be closed unless k = 0. Similarly, the two-point probability function V(9k, 
is also uniform unless k = 2k'. At first glance, this conclusion seems to contradict the reported 
non-uniform distribution of the phase difference of neighboring wavevectors D k = 9 k+ /\ k ~ @k m 
iV-body data (Scherrer, Melott, & Shandarin 1991; Coles & Chiang 2000; Chiang 2001; Chiang, 
Naselsky & Coles 2002; Watts, Coles, & Melott 2003), where Afc is a fixed small vector. The same 
arguments are also applied to higher-order approximations, so that the phase correlations between 
neighboring wavenumbers should not appear even in strongly non-Gaussian fields in a statistical 
sense. 

To resove this puzzle, it is useful to consider the conditional probability function given a Fourier 
coefficient of a small wavenumber a& k . I n the first-order approximation, the joint probability of 
having phases 9 k , #fc+Afc with fixed a Ak is given by 

V (0k,0k+Ak |c*Afc) oc 1 + ^|a A fe| cos {9 k+Ak -9 k - 9 Ak )p^(Ak,k), (12) 

which arise the non- uniform distribution pattern of the phase difference D k . The pattern depends 
on the fixed phase 9& k , which means the pattern varies from sample to sample, and this is exactly 
what is reported in the iV-body analyses. The pattern of the phase difference should be significant 
for red power spectrum, which is also consistent with iV-body analyses. The functional form of 
equation (12) also agrees with the iV-body analysis (Watts, Coles, & Melott 2003). The statistics 
of phase difference is thus the manifestation of the large-scale patterns of individual realizations. 
The position of the trough in the distribution of the phase difference corresponds to the phase of 
the mode Afc, and the degree of deviations from the uniform distribution depends on the specific 
amplitude of the mode Ak and also on the normalized bispectrum. 
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3. A numerical demonstration 

The equations (10) and (11) relate the bispectrum to the distribution of the phase sum 9 k +0 k i — 
k+k i. To see if this kind of phase information is practically robust, we numerically examine simple 
examples of non-Gaussian fields. Instead of examining cosmological simulations, the following 
simple example is enough to compare the numerical phase distributions and theoretical predictions. 
Series of non-Gaussian fields are simply generated by exponential mapping of a random Gaussian 
field: 

f(x) = exp (g<f>(x) - g 2 /2) - 1, (13) 

where is a random Gaussian field with zero mean, unit variance, and g is the non-Gaussian 
parameter. We simply take a flat power spectrum for the Gaussian field (p. The field / has zero 
mean and variance (f 2 ) = exp(g 2 ) — 1, and is called the lognormal field (Coles & Jones 1991). 
This field has quite similar statistical properties to gravitationally evolved non-Gaussian fields and 
approximately follows the hierarchical model of higher-order correlations. The parameter g controls 
the non-Gaussianity, and the random Gaussian field is recovered by taking the limit g — > 0. The 
random field / is generated on 64 3 grids in a rectangular box with the periodic boundary condition. 

In Fig. 1, the distribution of the phase sum 6 kl + 6 k2 — 9 kl+k , 2 is plotted for a binned con- 
figuration of the wavevectors, \k\\ = [0.4,0.5], [A^l = [0.5,0.6], #12 = [50°, 60°], as an example, 
where 6*12 is the angle between k\ and k2, and the magnitudes of the wavenumber are in units of 
the Nyquist wavenumber. The phase sum is averaged over the wavevectors in a configuration bin. 
The points represent the distributions of the phase sum in each realization. Poisson errorbars are 
smaller than the size of the points. The normalized bispectra p( 3 \k\, ^2) are numerically evaluated 
from each realization, which are used to draw the theoretical curves in the first-order approxima- 
tion of equation (11). There is not any fitting parameter at all. The agreement is remarkable 
in weakly non-Gaussian fields. When the non-Gaussianity becomes high, the data points deviate 
from the first-order approximation, and the distribution of the phase sum is sharply peaked at 
#fci + #fc 2 ~~ @k 1 +k 2 = mod 2-7T. Up to g ~ 3.0, or (f 2 ) 1 ^ 2 ~ 100, the distribution of the phase sum 
is accurately described by the first-order approximation, and is determined only by the normalized 
bispectrum. Even though the non-Gaussianity g ~ 3.0 on scales of the Nyquist wavenumber is 
beyond the perturbative regime, the normalized bispectrum on scales of the presently tested con- 
figuration is still within the perturbative regime, p^ ~ 0.25. This means that the phase sum is well 
approximated by first-order formula of the present work even when the field is strongly nonlinear 
in dynamics, as long as the parameter P{k)/V on the relevant scales is small. Increasing the power 
on relevant scales and/or decreasing the volume drive the phase correlation large, due to the fact 
that the phase correlations are particularly dependent on significant features in the sample. 
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ki + k2 - ki + k2 mod 2 77 

Fig. 1. — The distribution of the phase sum for a particular configuration of wavevectors. Five 
non-Gaussian fields are shown, where g is the non-Gaussian parameter and is the normalized 
bispectrum for the particular configuration. Theoretical predictions by normalized bispectra in the 
first-order approximation are shown by solid curves. 
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4. Summary 

The structure of the phase correlations in non-Gaussian fields is elucidated. The method to 
relate the joint distribution of phases to polyspectra is newly found and developed. The distribution 
of the phase sum of closed wavevectors is represented by the polyspectra. We found the statistics 
of the phase difference reflect the particular phase of the mode within an individual sample. The 
distribution of the phase sum of three or more modes carries the statistically useful information. 
The understanding of the phase correlations in non-Gaussian fields is now reached unprecedented 
level in this Letter, so that many investigations to make use of the phase information will be 
followed, such as the analysis of the non-Gaussianity of all kinds of cosmic fields, the nonlinear 
gravitational evolution of the density fields, the biasing and redshift-space distortion effects on the 
galaxy clustering, the primordial non-Gaussianity from inflationary models, and so forth. One may 
also hope that phase information can be useful in statistical analyses of all kinds of non-Gaussian 
fields, from various phenomena of pattern formations to human brain mapping, etc. 

I acknowledge support from grants MEXT 13740150. 
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